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able of contents 


§ #8 


D cL BRAT ons 
OTSSPOWCGCG_R3 = G COMPLEX*16 ** G COMPLEX®16 
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TS$POWCGCG = G COMPLEX*16 #* G COMPLEX®16 routine 16-SEP=1984 01: AX/VMS Macro_v04-00 p 
1802 routine TESEp=198 Wiovise EATMRTL. eeedorePcurcc.maRs1 2% Cty 


«TITLE OTSSPOWCGCG - G COMPLEX*16 ** G COMPLEX*16 routine 
«IDENT /1-002/ ; File: OTSPOWCGC.MAR Edit: SBL1002 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION,” MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 
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* THIS SOFTWARE is FURNISHED UNDER A LICENSE AND MAY BE USED AND gt ht 
* ONLY IN ACCORDANCE TERMS OF SUCH LICENSE AND WITH THE 
* INCLUSION OF THE ABOVE re NOTICE. THIS SOFTWARE OR ANY OTHER 
* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 

DIGITAL ASSUMES NO RESPONS . 


ITY FOR THE USE OR RELIABILITY OF ITS 


IBILITY 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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: ABSTRACT: 
: G COMPLEX*16 base to G COMPLEX*16 power giving G COMPLEX*16 result. 
: ENVIRONMENT: User Mode, AST Reentrant 


; AUTHOR: Steven B. Lionel, CREATION DATE: 20-July-1979 
; MODIFIED BY: 


; 1-001 - Original. Adapted from OTSS$POWCC version 1-003. SBL 20-Jul-1979 
; 1-002 - Use general mode addressing. SBL 30-Nov-1981 
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; FACILITY: Language support Library - user callable 
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- G COMPLEX*16 ** 
DECLARATIONS 


6 10 
G COMPLEX*16 routine . 


«SBTTL DECLARATIONS 
INCLUDE FILES: 


EXTERNAL DECLARATIONS: 


-DSABL GBL 

-EXTRN MTHSCGEXP 
-EXTRN MTHSCGLOG 
-EXTRN OTSSMULCG_R3 


MACROS: 


: EQUATED SYMBOLS: 


base = 4 
exp = 20 


t OWN STORAGE: 


> PSECT DECLARATIONS: 


1986 91:39:52 


AX/VMS Macro V04-00 


Complex exponentiation 
Complex logarithm 
Complex multiplication 


base input - by value 
exponent input - by value 


-PSECT _OTSSCODE PIC, USR, CON, REL, LCL. SHR, - 
EXE, RD, NOWRT, LONG 


r 
MTHRTL.SRCJOTSPOWCGC .MAR; 1 


Page 
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OTSSPOWCGCG = G COMPLEX*16 ** G COMPLEX*16 routine 16-SEP-1984 01:56:07 VAX/VMS Macro v04-00 
1800 OTSsPOUCEEG RS - G COMPLEX*16 ** G COMPL ~$eE 38 94 539:95 MTHRTL.SRCJOTSPOWCGC .MAR; 1 
8 aoe -SBTTL OTSSPOWCGCG_R3 = G COMPLEX*16 ** G COMPLEX*16 
p ; FUNCTIONAL DESCRIPTION: 
j : OTSSPOWCGCG_R3 evaluates the result of taking a complex base 
; to a complex power. The ANS FORTRAN X3.9-1978 standard defines 
¢ : complex exponentiation as: 
6 : ; x t* y = CEXP( y ® CLOG(x)) 
° 8 ; where x and y are type G COMPLEX*16. 
00 30 ; The arguments of OTSSPOWCGCG_R3 are CALL BY VALUE. 
9893 % : CALLING SEQUENCE: 
4 Be ; power.wgc.v = OTSSPOWCGCG_R3 (base.rgc.v, exponent.rgc.v) 
0000 3 : INPUT PARAMETERS: 
0000 98 ; Both base and exponent are G COMPLEX*16 numbers, each consisting 
0000 99 ; of a G REAL*8 real part and a G REAL*8 imaginary part. Both are 
0000 100; CALL BY VALUE. 
0090 101; 
0000 108 3; IMPLICIT INPUTS: 
0000 103; 
0000 104; NONE 
0000 105; 
0000 196 ; OUTPUT PARAMETERS: 
0000 107; 
0000 108 ; NONE 
0000 109; 
0000 110 ; IMPLICIT OUTPUTS: 
0000 111; 
0000 6 3 NONE 
0000 113; 
494 ie 3 FUNCTION VALUE: 
$600 ig : The G COMPLEX*16 (REAL*8, REAL*8) result of taking the 
0000 117; COMPLEX base to the COMPLEX exponent power is returned 
0000 118 3 in registers RO-R3. This is a violation of the VAX 
0000 119; cath ee standard, but is excused for compiled code 
0000 ! ? 3 support routines. 
$8 : > SIDE EFFECTS: 
099 ' : : Modifies registers RO-R3. 
3» : § : Possible error signals are: 
06 1 3 ; MTHS_INVARGMAT if base is (0.,0.). 
00 129 ; MTHS_FLOOVEMAT if floating overflow occurs. 
0 130 ; MTHS_SINCOSSIG if absolute value of the imaginary part of 
8 131; (exponent * CLOG(base)) > PIl#2**30. 
; ' § ; S$S$_ROPRAND if reserved floating operand is fetched. 
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OTSSPOWCGCG = G COMPLEX*16 ** G COMPLEX*16 routine -SEP=-1984 01:56:07 VAX/VMS Macro v04-00 Pa 4 

1-002 OTSSPOWC EG RS = G COMPLEX*16 ** G COMPL aie 94:39:83 MTHRTL.SRCJOTSPOWCGC .MAR; 1 - (4). 
188 «ENTRY OTSSPOWCGCG_R3, “M<> ; disable integer ovflo 
1 3 MTHSFLAG_ JACKET 3 establish math error handler 


6D  Q0000000°GF 9 MOVAB G*MTHSSJACKET_HND, (FP) 
set handler address to jacket 


handler 


SoOoOooooo 


COOQOQooOooooooooo 


L 
G 
0000 : 
9 
9 
83 1 
9 129 : Get complex logarithm of base 
= 7 & 3 146 7 SUBL2 #16, SP 3 return complex on stack 
04 AC oF OOC 14 PUSHAL base (AP) 3; address of base 
04 as F QOOF 144 PUSHAB ) 3; address of result 
00000000 ' GF 0 FB big 146 CALLS #2, G*MTHSCGLOG ; (SP) gets LOG(base) 
019 147 ;+4 
019 198 : CLOG(base) is at (SP). Multiply by exponent. 
BI3 133 3 Do multiplication out of Line. 
7E 1C AC 7D $019 151 ° MOVQ exp+B(AP), -(SP) ; Put exponent on stack 
7E 14 AC 7D 0010 136 MOVQ exp(AP), =(SP) ; CLOG(base) is already there! 
00000000'GF 08 FB 0021 15 CALLS #8, G*OTSSMULCG_R3 : RO-R3 gets CLOG(base) * exp 
Sh 1, 
0 138 : Now compute CEXP(product) 
$58 188 
ee ae a ee Jadigiennelsccenes 
5E 10 C2 QO2E 161 SUBL2 #16, SP ; Make room for result 
10 AE SF 0031 16¢ PUSHAB 16(SP) : Address of product 
04 AE OF te 16 PUSHAB 4(SP) 3; Address of result 
00000000 ' GF 02 FB 7 164 CALLS #2, G*MTHSCGEXP ; Result is at (SP) 
50 BE 7D OO3E 165 MOVa (SP)+, RO : Pop result into RO-R3 
52 8E 7D Bpeg 166 MOVQ (SP)+, R2 
Bae 168 RET : all done, exit 
045 169 END 
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- G COMPLEX*16 ** G COMPLEX*16 routine 


OTSSPOWCGCG 16-SEP-1 1: 
Symbol table 66-SEP-19 1984 oH 3$: 4 
BASE 4, 
exp = 900000 
be 9 yoy HND weeenere =X 1 
THSCGEXP rereerer x 
HSCGLOG eereerer x 
OTSSMULCG_R3 eeterere § Xx 
OTSSPOWCGEG_R3 00000000 RG 1 
¢oemmweeceenenwe ih 
: Psect synopsis ! 
Pees aewrescrenecan ceed 
PSECT name Allocation PSECT _No. Attributes 
- ABS. 00000000 ( g-} oO f 0.) NOPIC USR CON ABS LCL 
-OTSSCODE 0000004 ( 69.) 1 & Wes PIC USR CON REL LCL 
boewrtoe enema newer amnennemeom$ 
: Performance indicators ; 
Phase Page faults CPU Time Elapsed Time 
Initialization 30. 00:00:00.04 00:00:01.13. 
zonnens processing ot 0:00: 3-7 88:88:83 -0¢ 
Syabot table sort 9 0: 6: 0.01 90:00:00.01 
Pass 2 4 0:00:00.46 00:00:02.51 
Symbol table output ; 0:00:00.01 00:00:00.01 
Psect synopsis output 0: B08 «OF ee eS 
Cross-reference output 3 00:00:00.0 0:00:00.00 
Assembler run totals 28 00:00:07 .82 00:00:10.86 


et Limit was 750 
bytes (5 pages) of virtual memory were used to buffer the intermediate code. 
con were 10 pages of symbol table space at ecense to hold 7 non-local and 0 local symbols. 
29 source Lines were read in Pass 1, producing 11 object records in Pass 2. 
page of virtual memory was used to "define 1 macro. 


The working ages. 


Macro Library name 

“S255SDUAZ8:(SYSLIBJSTARLET.MLB:2 itt o- 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


NOSHR ae 
SHR 


AX/VMS 
MTHRTL. SRC 


cro v04-00 


Page 


€40 TSPOWCGC MAR; 1 


E NORD NOWRT NOVEC BYTE 
RD NOWRT NOVEC LONG 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWCGC/OBJ=O0BJ$:OTSPOWCGC MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MS 


; 
EE 
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